Lemaitre-Tolman-Bondi model of the universe 
and Hawking radiation of a dynamical horizon 



The paper deals with anisotropic spherically symmetric Lemaitre-Tolman- 
Bondi model of the universe bounded by the apparent horizon. Using Hamilton- 
Jacobi method for both massive and massless test particles, we are able 
to show that the temperature associated with the apparent horizon due 
to Hawking radiation is the usual Hawking Temperature. Also the usual 
Lemaitre-Tolman-Bondi model and the transformed r-gauge Lemaitre-Tolman- 
Bondi model are shown to be equivalent. Finally, the results of Tunneling 
approach agree with the Hamiltonian-Jacobi method. 

1 Introduction 

Hawking radiation(HR), a quantum description of a black hole(BH), is intimately related 
to the existence of an event horizon of a BH. The discovery of HR [Ij completes the cycle 
to describe BH as a thermodynamical object, i.e., BH behaves like a black body which 
emits thermal radiation with temperature proportional to its surface gravity and then 
Bekenstein [2: formulated the entropy proportional to the horizon area of the BH. This 
discovery also leads to view general relativity (GR) from different perspective. The inter- 
relation between gravity and thermodynamics opens new avenues in GR. The motivation 
comes from the equivalence of Einstein equations with the first law of thermodynamics 
{5Q = TdS), first illustrated by Jacobson ^ and subsequently, it is extended to non- 
Einstein gravity [4] O (6] [7] . At present, it is commonly believed that gravity might be 
originated from the thermodynamics of the unknown microstructure of space-time. Al- 
though, a self-consistent full quantum theory of gravity is not yet formulated but still 
BH thermodynamics is playing the role of a bridge to combine general relativity with 
quantum mechanics. 

The derivation of Hawking that "BH evaporates particles" [TJ [S] was based on the 
quantum field theory. Subsequently, Hartle and Hawking [9] derived the BH temperature 
semiclassically, by using Feynmann path integral. Due to mathematical complexity of the 
above procedures, semiclassical approaches [101 1111 1121 113) were developed for studying 
BH radiation. Basically, these semiclassical techniques were classified into two approaches 
— one, due to Parikh and Wilczek [TOl[TT] and the other developed by Padmanabhan et. 
al. [HIIS]. 
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In the first approach, the main ingredient is the consideration of energy conversion in 
tunnehng of a thin shell from the hole and is known as the radial null geodesic method. 
Here the imaginary part of the action for the s-wave emission is related to the Boltzmann 
factor for emission to relate with Hawking temperature(HT). On the other hand, the 
second method is known as complex paths method. Here the action S{r, t) for a single 
scalar particle is obtained by solving the Klein-Gordon(KG) equation in gravitational 
back ground using the Hamilton-Jacobi(HJ) approach. Then HT is obtained using the 
"The principle of detailed balance" [HI [T5] . 

So far, most studies of HR effect deal with static space times, where there exist an 
event horizon. The global concept of event horizon is used to define the HT. But in 
dynamical space time, difficulty arises due to non existence of event horizon even locally. 
However, using the tunneling approach of Parikh-Wilczek [10], recently Hayward et. al. 
[14j proposed a locally defined HT for dynamical BH. Subsequently, Cai et. al. [15] 
have shown HR from the locally defined apparent horizon of the FRW universe where 
HT is measured by an observer using the Kodama vector |16l 1171 ITSj inside the horizon. 
Thus, at present one of the important issue in GR is to formulate the thermodynamics 
of the dynamical space times and its relation with gravity. Formulation of HR is one 
step forward towards this goal. In the present work, we deal with spherically symmetric 
LTB model of the space time and derive the HT using H-J method both for massless and 
massive particles. The semiclassical tunneling approach has been described and it results 
identical HT on the horizon. Finally, the paper ends with a brief discussion. 

2 Lemaitre-Tolman-Bondi model and Hamilton- 
Jacobi method for massless particle 

The inhomogeneous spherically symmetric LTB space-time model is described by the 
metric ansatz in a co-moving frame as 

ds^ = -dt"^ + Trrdr^ + {dO"^ + sin^ ed^)^) (1) 

l + /(r-J ^ ' 

where R = R{r,t) is the area radius of the spherical surfaces and /(r)(> —1) is the 
curvature scalar classifies the space-time as 

(i) bounded, if, -1< f{r) < 

(ii) marginally bounded, if, /(r) = 

(iii) unbounded if /(r) > 0. 

The Einstein field equations for the space time model can be written as |19ll20ll2ni22| . 



and the evolution equation for R is 

2RR + if + 8nGpR'^ ^ f{r). (3) 

Here the mass function 

F{r,t)=R{R'~f{r)), (4) 

is related to the mass contained within the co-moving radius V. 

The universe is assumed to fill with perfect fluid with energy-momentum tensor 

T^f — {p + p) Ufj,u„ + pg^i, (5) 
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where the fluid 4- velocity u'^ has normahzation u^u^ = — 1 and p, p are respectively 
the usual matter density and pressure of the fluid. The energy momentum conservation 
relation T.^"^ — demands the thermodynamic pressure to be homogeneous (i.e., p' = 0) 
and the evolution of matter density is given by 

p + 3H{p + p) = (6) 

where, H = ^ ^-S^ + 2-^^ is the Hubble parameter, 'dot' and 'dash' over any quantity 

stand for differentiation w.r.t. 't' and 'r' respectively. Now splitting the metric ansatz 
^ on the surface of the 2-sphere and on the 2-D hyper surface normal to the 2-sphere as 

ds^ = habdx'^dx^ + R^d^ (7) 

the dynamical apparent horizon is characterized by |23l 124) 



h^^daRdiR = (8) 

where hah = diag i+f{r) ) ^-E) metric normal to the 2-sphere. Thus the 

spherical surface of radius R = Ra corresponding to the apparent horizon (i.e., marginally 
trapped surface) satisfies 

Ra = F{r, t) and R\^1 + f{r) (9) 

On the other hand , trapping horizon {Rt), a hyper surface foliated by marginal spheres, 
is characterized by [231 124j 



d+Rr = , i.e., Rr = a/1 + f{r) (10) 

where, 

a. = -^(a.T^i5^9.). (11) 

are the null vectors normal to the 2-sphere. 

Hence both the horizons coincide for LTB model as in FRW space-time. 

Now, the Misner-Sharp gravitational mass (in units G = 1) is defined as [241 125) 

m(r,t) = ^ (l - h^'daRdbR) (12) 

one may note that this mass 'm' is an invariant quantity on the 2-D hyper surface nor- 
mal to the 2-sphere and m = niH = on the horizon. Further, one can introduce 
another invariant scalar associated with the horizon in the normal hyper surface, known 
as dynamic surface gravity and is defined as [H] 

KD = —^da h'^^'dbR) \h (13) 



which for the given model has the expression 



^ [2^ {"^' (^^') + \^'^'^] 



(14) 
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For the LTB model with decomposed metric (O, the Kodama vector k has components 



n''{r,t) = 



-.e'^^'dbR 



K = = 



(15) 



R 



So = T?'^ — 1 — /(y), i.e., Kodama vector is time-hke, null or space like for inside, 

on the surface or outside the apparent (i.e. trapping) horizon respectively. Note that 
Kodama vector is very similar to the time like Killing vector for stationary BH space 
time and consequently, an invariant energy associated with a particle is defined by the 
scalar on the normal space as [26j 



U) = —K^dal ■ 

Here the classical action / of the massless particle satisfies the HJ equation 



h'^'^daldtl = 



(16) 
(17) 



(diy 

Hence solving (|16p and (I17p we have 



1 + I (dl 



dr 



ujR' 



and ^ 
dt 



{7?-^l + /(r)} 



(18) 



where there is a pole at the horizon. 

The full classical action of an outgoing massless particle is 

I — daldx"^ 



(19) 



where 7 is an oriented curve with positive orientation along the increasing values of 
x°' = {t,r). As for massless particles the radial motion is along a null direction so from 
the metric we have 

p/2 

= ds^ = -dt^ + rr—^rr^dr'^ 



1 + f{r) 



R' 



dt = ±^^^dr, (20) 

for outgoing/ingoing particles respectively. Hence using equation p8[) and 1)20^ in the 
equation p9|l we have for outgoing particles 



1 = 2 dr {dri) = -2 



LUR' 



dr 



n + f(r)}(, R \ 

I 0T7m J 



(21) 



Now expanding G{r,t) = 1 
direction we obtain 1271 



, ^ in the neighborhood of the horizon along a null 



Gir,t) 



R 

Vr+7' 



-.At- 



R' 



Rf'(r) 



Vl + .Z 2(l + /)5 



Ar 



+ 



4 



RR' 



R' 



Rf'ir) 



(1 + /) VTTJ 2(1 + /)^ 

27?' 



Aria + . 



(1 + /) 



Then from ((2T} 



dr- 



(22) 



(23) 



{r — th — iO) ' 

which has a simple pole aX r = rn- Using Feynmann's ie-prescription, the imaginary part 
of the action (the real part has no physical consequence) can be written as [27] 

-KbJH 



Iml ■ 



(24) 



Hence, one may interpret T — > as the dynamical temperature associated with 

LTB space-times. 

Moreover, we consider LTB tunneling computation in the co-ordinate system (f , t, 6, <j) 
where f = R. The metric ansatz (in r— gauge) becomes 



ds^ = -Adt'^ - 2Bdfdt + C{df f + r^dnl 
= Qabdx'^dx'' + r^dQ,2 



(25) 



where, 



4 = 1- 

B = 
and C - 

the horizon is located at x = where 



R^ 



1 + /W 

R 



1 + f{r) 
1 



1 + f{r) 



= A(f,t) 
B{r,t) 
= C(f,t) 



ab i~i ~. f~i ~ rr 

X = g OaTObr = q = 



A 



AC + B2 



i.e.. Ah = gives the horizon. 

The Kodama vector is given by 



1 



^VB^+AC 
and the invariant energy has the expression, 

dtl 



VB^ + AC 

The dynamical surface gravity on the horizon is 



1 

2B3 



B+-AC 



(26) 



(27) 



(28) 



(29) 



The HJ equation for a massless particle along a radial trajectory has the explicit form 



C {dtif - 2B (dtl) {dfl) + A {dflf = 



(30) 
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Now integration over the temporal coordinate gives real contribution to the particle ac- 
tion (which has no physical significance), an imaginary contribution comes only from 
integration along the radial direction . So, eliminating dtl between equation (|28p and 
H30|l we obtain 

dfl = ^ '-^Jb-^^AC (31) 

For the metric (125^ . by making a null expansion on the horizon we obtain 

= ds^ = -2dtdf + 



and 



\ 2RRdt I 2RR' f{r)R^ 



l + f{r)J l + f{r) \l + f{r) (l + f{r)) 
\ RR 2RR' f'{r)R^ 1 



\2 



df 



Hence 



= 2KDH{f-fH) (32) 

/../^-/.P^ii^±^^!l^^-^ (33) 



Thus we have the same result. 



3 Massive particle 



We consider a particle of mass 'm' moving radially in the space time described by the 
metric ([25}. The HJ equation is 

g^^d^Sd^S + = (34) 
where the action S can be written as 

5= / dt+ Kfdf (35) 



Here, Kf, the radial momentum associated with the test particle and from equation 

as 
at 



^ is related to the invariant energy oj as 



^ = -lj^/B^ + AC (36) 

hence we write, 

S = - j ^VB'^ +ACdt + J Kfdf (37) 



Now, the explicit form of the HJ equation (I34|l for the action (I37|) and the metric (I25|) is 
given by 
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AKf - 2ujKrB^/B^ + AC + (S^ + AC) (m^ - Cuj^) = 

which gives 

VWTac 



^ujB ± ^Lj2 (B2 + 4(7) - ^m^l (38) 



As the observer is inside the apparent horizon so considering the incoming mode (negative 
sign in the expression for Kr; positive sign corresponds to an outgoing mode), the imagi- 
nary part of the action is 

ImS = J + ^'^ |cjB ± y/io^ (B2 + AC)- Am^l df (39) 

It is to be noted that the integrand has a pole at the horizon, so performing the contour 
integration we have [IF| 

ImS = Trrffo; (40) 
In tunneling approach, using WKB approximation the emission rate F (for particles tunnel 
from outside to inside the horizon) is the square of the tunneling amplitude TT5' 

i.e., T (X exp{—2ImS} — exp{—2nrHijj} (41) 

which can be casted in the form of a thermal spectrum F ~ exp { — ^ } provided T = ^^l^ , 
the HT. 

Thus an observer, inside the apparent horizon will see a thermal spectrum with HT 
when particles tunnel from outside to inside the horizon. 



4 Tunneling approach 



The basic idea in tunneling method is that following the standard approach, i.e., in the 
semiclassical approximation (i.e., WKB approximation), the emission rate of tunneling of 
a massless particle across the horizon can be related to the imaginary part of the action of 
the system. In the s-wave approximation, the particles are considered as massless shells, 
moving along a radial null geodesic. 

So for the metric (|25|l the radial null geodesic is characterized by 



.^B±VWTAC (^2) 
o 

where as before '+' or ' sign indicates an outgoing or incoming null geodesic. Due to the 
tunneling of the particles from the outside to the inside the horizon, here we consider only 
an incoming geodesic. As we have seen that in the present context we need the imaginary 
part of the action produced by the tunneling particles (remaining part is always real) i.e., 
particles tunneling through a barrier (the classically forbidden region) so we obtain [T5] 

ImS = Im / pf-dr — Im I / dp'fdf (43) 
Jfi J f-a Jo 

Here, the particle with radial momentum pf tunnels from the initial position just 
outside the horizon to the final point at f/ which is a classical turning point, i.e., in the 
semiclassical analysis, the trajectory can represent a classically allowed motion 



dH 

dp,-- 



(44) 
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where the Hamiltonian H is the generator of the cosmic time 't' . Then using equation 
(flil) in equation (gSJ we write [TS] 

ImS = Im f df f = Im f ^lo\/ + AC = —uilm f d fC\/B + _ ^f.^ 

Jf, J f Jr, r J,^ {VB^ + AC~B} 

(45) 

One may note that to perform the integration over the Hamiltonian H we get the energy 
as ujy/ + AC, as measured by an observer with the Kodama vector. Then using the 
emission rate (|4ip we have 

which is the HT as derived in the previous section. 



5 Discussion 



The paper deals with thermodynamics of the inhomogeneous LTB model of space-time. 
In the HJ formalism, we have considered the tunneling of both massless and massive 
particles, and the HT is measured by an observer with the Kodama vector (eauation (|15|l ') 
inside the apparent horizon (the trapping horizon). Thus apparent horizon of LTB uni- 
verse is related to the HT in contrast to the BH case when HR is associated with the event 
horizon. Also the present work supports the results of FRW universe p3] in the literature. 
Also making a co-ordinate transformation (given by the metric equation (|25p l we reach 
to the same conclusion, indicating HR independent of coordinate choice. Subsequently, 
tunneling approach (proposed by Parikh and Wilczek) has been described for massless 
particle to show the equivalence between the two semiclassical approaches. Finally, it will 
be interesting to improve our results by taking into account of back reaction of HR as in 
case of BH, consideration of back reaction results in a non-thermal spectrum of HR and 
as a result the issue of information loss may be overcome. Also for future work, it will 
be interesting to consider higher order correction terms and examine whether quantum 
correction to Hawking temperature [28] may be obtained. 
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